Large Scale Intermittency in the Atmospheric Boundary Layer 
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We find actual evidence, relying upon vorticity time series taken in a high Reynolds number 
atmospheric experiment, that to a very good approximation the surface boundary layer flow may be 
described, in a statistical sense and under certain regimes, as an advected ensemble of homogeneous 
turbulent systems, characterized by a lognormal distribution of fluctuating intensities. Our analysis 
suggests that usual direct numerical simulations of homogeneous and isotropic turbulence, performed 
at moderate Reynolds numbers, may play an important role in the study of turbulent boundary layer 
flows, if supplemented with appropriate statistical information concerned with the structure of large 
scale fluctuations. 
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A major difficulty in dealing with boundary layer flows 
at high Reynolds numbers is that they cannot be straight- 
forwardly modeled within the theory of homogeneous 
isotropic turbulence, rendering unlikely, in principle, an 
application of the general results of the latter in the con- 
text of the former. In fact, essentially all the symmetry 
properties of the evolution equations break down close to 
the boundaries due to the intermittent production of a 
whole "zoo" of flow phenomena, like low speed streaks 
and a number of vortex structures. The lack of homo- 
geneity and isotropy is anagously observed far enough 
from the walls, where the transition to the outer laminar 
flow takes place, a region of strong entrainment and high 
intermittency factor, as discussed several decades ago by 
Klebanoff in his benchmark work (and subsequent pa- 
pers) on turbulent boundary layers [l| . 

It has been long hypothesized, however, that in a typ- 
ical turbulent boundary layer problem, as in the flow 
over a flat or rough surface, there is an intermediate 
range of normal distances from the boundary - the log- 
arithmic layer - where the fundamental symmetries of 
the Navier-Stokes equation are approximately restored 
at small scales, yielding a stage for tests of the statistical 
theory of turbulence. A recent investigation of related is- 
sues is provided by Sreenivasan et al 2\ , from the analysis 
of the time series produced by hot-wire anemometry in 
a stable atmospheric boundary layer. It is worth noting 
that in atmospheric experiments it is usual to get samples 
where the flow is at best approximately statistically sta- 
tionary, which makes the connection with the physics of 
homogeneous and isotropic turbulence not obvious at all. 
Actually, the experience shows that while the original ve- 
locity signal can be used to check the Kolmogorov's four- 
fifths law or the scaling exponents of structure functions, 
for instance, the lack of stationarity has to be carefully 
accounted for in the study of probability distributions of 



local fields. As a pragmatical solution of such a "large 
scale intermittency" problem, the authors of Ref. Q have 
retained from the rough anemometric data only the sam- 
ples that would lead to statistically stationary regimes. 
It turned out, a posteriori, that their prescription worked 
consistently well. 

The fact that several scaling features of homogeneous 
turbulence are observed without any further handling of 
the velocity signal seems to indicate that the underlying 
boundary layer flow could be modeled, in a first approx- 
imation which disregards shear effects, as an ensemble of 
homogeneous turbulent systems, collectively advected by 
the mean local velocity U . Each element of the ensemble 
would correspond to a flow with a definite value of the 
turbulent intensity / = u rms /U . Thus, if we are inter- 
ested to model fluctuations of some observable defined at 
length scale £, 0(£,t), we may write 



0(t,t) = x(t)d(£,t) 



(1) 



where 0(£,t) denotes the observable fluctuation associ- 
ated to an arbitrary homogeneous and isotropic turbulent 
flow, and x(t) is an independent random function of time, 
which accounts for the fluctuations of the rms values of 
0(£, t). In effect, x(t) may be thought to play the role of 
a positive enveloping function which modulates the faster 
fluctuations of 0(£ , t). 

As an illustration of the kind of modeling we have in 
mind, take the case of longitudinal velocity differences 
ve(t) = vi(r + ix, t) — vi(r, t). We have, from |T|), 



v e (t) = x(t)vi(t) . 
It is clear, therefore, that 

(vj) = S q (i) = {x*){& e ) 



(2) 



(3) 



that is, the structure functions S q (£) and S q (£) depend 
on the length scale £ exactly in the same way. On the 
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other hand, it is not difficult to find that the probability 
distribution functions (pdfs) of V£ and vi are, in general, 
completely different: we have 



P{vt) 



dx\x\ 1 p(v e /x)f(x) 



(4) 



where f(x) is the pdf for the random variable x intro- 
duced in |T]), and p(-) and p(-) refer to the pdfs of velocity 
differences in boundary layer and homogeneous isotropic 
flows, respectively. 




'////////////////////////////////////////////////////////////////////////////. 

WALL 



FIG. 1: Coherent structures generated close to the wall are 
transported to the bulk of the flow. These configurations are 
subject to random advection and instabilities causing them to 
disrupt, grow, and burst as the evolution proceeds. Approxi- 
mately homogeneous and isotropic turbulence regions of size 
of order L are produced with random energy transfer rates e. 

Our central aim in this note is to discuss the statis- 
tics of the random enveloping function x(t), through the 
analysis of the vorticity time series obtained from an at- 
mospheric experiment carried out at very high Reynolds 
number Q ■ The measurements were performed with a 20 
hot-wire probe, which incorporated specific design fea- 
tures appropriate for the particularities of the field ex- 
periment. The data was collected at a sampling rate 
of 10 KHz (which was high enough to resolve the dissi- 
pative Kolmogorov scale) in a tower 10 m high, placed 
over a grass-covered flat surface. The total time length 
of the velocity signal is 15 minutes, corresponding, by 
the mean wind velocity, to 6.3 Km. The anemometer 
set was calibrated at the measurement position, to avoid 
possible perturbations caused by its manipulation and 
the variation of environmental factors. The Taylor-based 
Reynolds number of the flow, observed in approximately 
neutral and stable conditions, is R\ ~ 10 4 ; for a more de- 
tailed account of the experimental and phenomenological 
parameters, see 0. The experiment outcome, then, con- 
sists of time series for the three velocity components and 
all the nine components of the velocity-gradient tensor 
as well. 

We intend to check, on the basis of purely heuristic 
and phenomenological arguments, whether the random 



variable considered in Eq. (QJ, is lognormally dis- 

tributed. A turbulent blob of size L and rate of energy 
dissipation e is produced in the boundary layer along a 
complex cascade directed towards larger scales, as shown 
in Fig.l. As a working hypothesis, an analogy with 
the multiplicative cascade arguments of the K62 phe- 
nomenology [illH can be drawn here, assigning lognormal 
fluctuations to e. Its random behavior would be the re- 
sult of sucessive disruptions of coherent vortex structures 
generated at the boundary, followed by their straining 
and transport to upper positions in the flow. Since x(t) 
may be essentially identified to fluctuations of the rms 
velocity, the usual expression 



(cL)t 



(5) 



taken from the statistical theory of turbulence [6J| , estab- 
lishes a connection between x(t) and e. Therefore, we 
expect that u rms , or, equivalently, x(i), will be lognor- 
mally distributed, at least approximately. 

Manifestations of lognormal statistics in the physics of 
boundary layer flows are not unusual. It is worth men- 
tioning that such distributions are found in turbulent 
boundary layers for: (i) the time intervals between ve- 
locity "bursts" 0|; (ii) the spanwise separation between 
streaks intermittently produced at the wall Q, and (iii) 
the small scale fluctuations of the dissipation field . An 
interesting question is if all of these instances of lognor- 
mality can be related to each other within the framework 
of some unifying description. 

Assuming that the ^-variable in Eq. |T]) is lognormally 
distributed is just half of the whole story. We also have 
to model the fluctuations of the small-scale observable of 
specific interest, which we take to be the vorticity field, 
k>i = CijkdjVk- There are experimental indications that 
even at moderate Reynolds numbers the enstrophy pdf 
(scaled to have unity variance) has, for a large range of 
enstrophy values, a turbulent asymptotic profile [lj| • Re- 
lying on isotropy, that observation suggests that a similar 
asymptotic behavior should hold for the vorticity fluc- 
tuations at Reynolds numbers which are not necessarily 
very high. In particular, it is likely that the vorticity pdfs 
obtained from direct numerical simulations, as the ones 



performed by Vincent and Meneguzzi ll| have asymp- 
totic shapes (we mean simulations with R\ ~ 150). 

We have defined, motivated by the above considera- 
tions, a smooth polynomial interpolation of the Vincent- 
Meneguzzi numerical vorticity pdf and used that interpo- 
lation to generate a stochastic process for <D, an arbitrary 
component of vorticity, by means of a standard Monte- 
Carlo procedure [l2| • More precisely, let i — 1,2... be an 
integer index and x(i) and ui(i) be independent random 
variables distributed according to the lognormal and the 
Vincent-Meneguzzi vorticity pdf, respectively, both with 
fixed variances. A stochastic process which simulates the 
vorticity fluctuations measured in the boundary layer is 
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given by the series defined by 



Lu(i) = x{i)Cj(i) . 



(6) 



The lognormality of the x-variable is implemented with 
the help of the mapping x = exp(y), where y is normally 



distributed with variance a y . Observe that 



H(n) = 



_ (uj 2n ) _ (x 2n ) {Co 2n ) 



{uj 2 ) n (X 2 ) n (0J 2 ) n 



_ _2n(n-l) 



■vH(n), (7) 



a result that gives a hint on why the hyperflatness fac- 
tors H(n) are reasonable higher in boundary layers than 
the ones typically found in homogeneous and isotropic 
turbulent flows. We have considered a stochastic process 
with 31 x 10 6 elements, for various values of a y . The 
choice o y = 0.51 leads to an excellent agreement with 
the empirical vorticity pdfs, as shown in Fig. 2. 
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FIG. 2: Vorticity histograms. The black line is the result 
of the simulated stochastic process; the underlying colored 
lines are associated to the components of the vorticity vector, 
measured in the atmospheric boundary layer. 

An important point concerns with the precision as- 
signed to a y , and how it affects the vorticity distribution. 
We have found that alternative profiles for the pdf of u> 
and appropriate redefinitions of a y lead to reasonable fit- 
tings as well. If, for example, uj is taken to be normally- 
distributed (which is obviously wrong) we get a y = 0.7. 
On the other hand, if we take a Student's t-distribution, 
which models reasonably well the vorticity pdf tails, as 
advanced in Ref. fl3j ]. a very good agreement is found 
again, this time with a y = 0.61. The uncertainty in the 
value of a y is partially due to properties of the lognormal 
distribution and partially related to the structure of Eq. 



@. We will not discuss these mathematical aspects in 
detail, but it suffices to note that a random variable x, 
lognormally distributed, can be decomposed arbitrarily 
as the product of two other independent random vari- 
ables, x\ and X2, both following lognormal distributions. 
This means that we may write 



x(i) = Xi(i)x2(i) 



to get, from ©, 



uj{i) = xi{i)[x 2 (i)Cb(i)} = xi(i)w(i) 



(8) 



(9) 



where x\(i) takes the place of x(i) and = X2{i)Co{i) 
takes the place of Cj{i) in Eq. ©. We find, thus, that 
the mere fitting of the vorticity pdf through the use of 
the stochastic process ^ has to be interpreted with care, 
even if the fitting is incredibly accurate. 

We conclude this note emphasizing that the problem 
of modeling fluctuations of vorticity - or any other small- 
scale observable - in the boundary layer is not so straight- 
forward as it could seem at first sight. While we have 
demonstrated that the use of the lognormal distributions 
and statistical modeling based on DNS lead to excellent 
results, more elaborate statistical information is in or- 
der. A promising approach is to find instances where the 
same enveloping function x(t), introduced in Eq. (TTJ, 
would modulate intermittent fluctuations of alternative 
observables. We suggest filtered velocity fields as good 
candidates. Nevertheless the usual gaussian behavior of 
velocity, it is known that a frequency-filtering procedure 
applied to the velocity time series produces non-gaussian 
pdfs [2| , which, hopefully, could be modeled with the help 
of the lognormal random functions advanced here. Since 
this an issue of far more detailed analysis, a comprehen- 
sive discussion is deserved for future work. 

During the completion of this paper, we have become 
aware of an article by F. Bottcher et al. [14j |. which has 
some overlap with the work presented here, mainly in 
the interpretation of boundary layer flows in terms of 
homogeneous isotropic turbulent ensembles. 

We thank Marie Farge for sending the numerical vor- 
ticity pdf files computed from the Vincent-Meneguzzi 
DNS. L.M. thanks partial financial support by CNPq and 
FAPERJ. 
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